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Abstract — Monitoring environmental changes from space re- 
quires extremely well-calibrated observations to achieve the nec- 
essary high accuracy and stability. The calibration differences 
between the Moderate Resolution Imaging Spectroradiometer 
(MODIS) and the Advanced Very High Resolution Radiometer 
(AVHRR) thermal bands provide a valuable quality assessment of 
the instrument performance. This letter compares the calibration 
differences between the Aqua MODIS and NOAA-18 AVHRR 
bands at 11.0 and 12.0 /xm using simultaneous nadir overpass 
observations obtained in nearly parallel orbits. Impacts due to the 
relative spectral-response differences between the two sensors are 
estimated by MODTRAN simulations with real-time atmospheric 
profiles of temperature, water vapor, atmospheric pressure and 
ozone, and surface skin temperatures. Results show that the tem- 
perature difference after the removal of atmospheric impacts is 
within 030 K (or 0.40% in radiance) across the effective calibra- 
tion range for the 11.0-/xm band/ channel. For the 12.0 - fim band, 
the differences are 0.40 K (or 0.50%) at the typical radiance and 
up to 0.70 K (or 0.90%) close to the maximum radiance, indicating 
an excellent calibration consistency between MODIS and AVHRR 
for both bands. 

Index Terms — Advanced Very High Resolution Radiometer 
(AVHRR), calibration, MODerate Resolution Imaging Spectrora- 
diometer (MODIS), MODerate resolution atmospheric TRANs- 
mission (MODTRAN), simultaneous nadir overpass (SNO), 
temperature. 


I. Introduction 

T HE Moderate Resolution Imaging Spectroradiometer 
(MODIS) is a key instrument of the National Aeronautics 
and Space Administration’s Earth Observing System (EOS) 
mission currently operating on both Terra and Aqua satellite 
platforms (launched in December 1999 and May 2002, respec- 
tively) in a near Sun-synchronous polar orbit [l]-[3]. MODIS 
is designed from a number of heritage sensors, including the 
Advanced Very High Resolution Radiometer (AVHRR), which 
has been in operation for nearly three decades on board the 
National Oceanic and Atmospheric Administration (NOAA) 
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polar-orbiting satellites [4], MODIS and AVHRR are the two 
most important sensors that make global observations to derive 
various products which monitor land, ocean, and atmospheric 
processes. Monitoring these environmental changes from space 
requires extremely well-calibrated observations to achieve the 
required high accuracy and stability [5]. 

The MODIS instrument has 36 spectral bands covering a 
wavelength range from 0.4 to 14.4 |Um with Earth scene ob- 
servations sampled at three nadir resolutions of 0.25, 0.50, 
and 1.0 km. The 16 thermal emissive bands (TEBs, i.e., band 
20-25 and 27-36) are calibrated on-orbit using a large aperture 
V-grooved blackbody (BB) controlled at a constant temperature 
of 290 K for Terra and 285 K for Aqua. The three AVHRR 
thermal infrared channels (3.7, 11.0, and 12.0 of channels 
3b, 4, and 5, respectively) are calibrated in-flight with a full- 
aperture honeycomb BB as an internal calibration target (ICT). 
However, the ICT temperature varies with instrument ambient 
environment, with temperature floating at around 290 K [6], 
[7]. The absolute radiometric calibration of the MODIS TEB 
and AVHRR infrared channels is related to the standards of the 
National Institute of Standards and Technology. 

Although the calibration of MODIS and AVHRR sensors 
each relies on measurements of the on-board BB and the 
cold-space background, the calibration methodology operates 
differently, depending on the design of the instrument. The 
independent calibration of thermal bands between MODIS and 
AVHRR provides a valuable quality assessment of the perfor- 
mance of each instrument. A direct comparison of tempera- 
ture obtained from near-simultaneous nadir overpasses (SNOs) 
shows that the differences are generally within 0.25 K at typical 
scene temperatures for the atmospheric window bands [8]-[10]. 
A clear dependence of the temperature difference on scene 
temperature has been found between MODIS and AVHRR 
and for most pairs of NOAA series AVHRR [11]. Since the 
MODIS spectral-band passes are significantly narrower than 
those of the AVHRR channels, temperature differences between 
the two sensors could be scene temperature dependent and 
atmospheric water vapor sensitive. However, the results of pre- 
vious MODIS and AVHRR comparison did not account for the 
impacts of atmospheric water vapor content. To better quantify 
the sensor differences caused by calibration, it is necessary 
to remove atmospheric associated impacts on the calculation 
of the temperature differences. In the past, the lack of real- 
time observations of atmospheric temperature and humidity 
profiles that coincided in time with satellite observations made 
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Fig. 1 . RSR of the two spectrally matched bands/channels (1 1 .0 and 12.0 /zm) of Aqua MODIS (3 1 and 32) and NOAA-1 8 AVHRR (channels 4 and 5). 


it difficult to evaluate such impacts since the atmospheric 
temperature and water vapor profiles change constantly with 
time and geolocation. Recent developments and improvements 
in the sounding product have allowed sensor-retrieved real- 
time surface temperatures and atmospheric profiles to become 
available. 

This letter uses observed temperature differences obtained 
from Aqua MODIS and NOAA-1 8 AVHRR SNO observa- 
tions to examine their calibration differences for thermal 
bands at 11.0 and 12 ^m. Impacts of atmospheric wa- 
ter vapor on the calculation of the temperature differences 
are simulated by MODerate resolution atmospheric TRANs- 
mission (MODTRAN4), Revision 1 [12], with real-time at- 
mospheric profiles obtained from the Atmospheric Infrared 
Sounder (AIRS) on the Aqua spacecraft. Currently, the lat- 
est version (MODTRAN5) is available (http://www.ontar.com/ 
Software/ProductDetails.aspx?item— modtran), and significant 
improvements are made for calculations of solar spectral ra- 
diances. Thus, it still justifies our use of MODTRAN4 in this 
letter for the atmospheric window spectral region. We selected 
Aqua MODIS and NOAA-1 8 AVHRR because their SNO 
events occur periodically in overlapping orbits with crossover 
angles at nearly 0°, thus providing extra opportunities to ex- 
amine and compare calibration stability of the two sensors. 
To ensure the quality of the comparison, only data collected 
over relatively homogeneous surfaces are used. Orbits over 
different latitude sections are selected to cover a wide range 
of scene temperature. The objectives of this letter are to: 1) 
provide MODTRAN simulations of MODIS and AVHRR real- 
time temperature differences caused by their relative spectral 
response (RSR) differences and 2) provide quantitative as- 
sessments of MODIS and AVHRR calibration differences by 
combination of SNO observations with the real-time simulation 
of atmospheric influences. 

II. Methodology 

The comparison between Aqua MODIS and NOAA-1 8 
AVHRR is conducted for two spectrally matched bands (Fig. 1), 
which are MODIS band 31 (11.03 ^urn) to AVHRR channel 
4 (10.79 fjm) and MODIS band 32 (12.02 /xm) to AVHRR 
channel 5 (11.92 /nn). The effective temperature (T), at which 
the underlying surface radiates as a BB, is used to evaluate the 
differences between the two sensors. To accurately compute T, 


the measured brightness temperature ( Te ) is calculated based 
on the central wavelength (A c ) of a given band 


T e 


Cj/K. 

ln[l + Cl /(LA»)] 


( 1 ) 


where ci and C 2 are the Planck equation constants and L is the 
radiance of a band. To determine T, a linear regression between 
T and Te is used 


T=(T e - A)/B. (2) 

Coefficients A and B for NOAA series AVHRR are pro- 
vided in the NOAA-KLM User’s Guide (https://www.ncdc. 
noaa.gov/oa/pod-guide/ncdc/docs/intro.htm), while coefficients 
for Terra and Aqua MODIS are derived using the same proce- 
dure [8], [9], [13]. Basically, coefficients A and B are deter- 
mined from the linear regression between Te and T, with T 
changing from 180 K to 320 K at an interval of 0.1 K. The 
theoretical value of Te is computed by convolving the RSR 
with the BB spectral radiance at T. The relative difference in T 
between MODIS and AVHRR is calculated by 

AT = TmoDIS — ^avhrr (3) 

where Tmodis and Tavhrr are the effective temperatures for 
MODIS and AVHRR, respectively. 

Theoretical values of AT due to the MODIS and AVHRR 
RSR differences are estimated from MODTRAN simulations. 
Our simulation uses AIRS-observed real-time temperature, 
pressure, humidity and ozone profiles, and skin temperatures 
obtained during the MODIS and AVHRR SNO periods. Thus, 
the atmospheric profile data are measured at the same time 
as the Aqua radiance measurements, while time lags in SNO 
data between MODIS and AVHRR are restricted to be within 
30 s. Each simulation provides a real-time top-of-atmosphere 
(TOA) radiance spectrum [L(A)]. Given the RSR for MODIS 
and AVHRR, a simulated band spectral radiance (L simu ) is 
determined using 

iSimu Sxi RSR(X)L(X) 

J^RSR(X). 

By replacing L in (1) with L simu and following the same 
calculation path from (1) — (3), the simulated relative difference 
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Fig. 2. Histogram of sampled 50-km areas along latitude. 

in T between MODIS and AVHRR can be determined. The 
remaining differences between the sensors observed and sim- 
ulated (i.e., theoretical values) AT are considered to be caused 
by relative calibration differences, i.e., 

AT'calib __ ^rpobs f \rp simu 


III. Results 

A. Observed and Simulated AT 

The observational data are collected from 13 Aqua MODIS 
and NOAA-18 AVHRR SNO events that occurred between 
2006 and 2010 in nearly overlapping orbits, meaning that the 
crossover angles between the two sensors are close to 0°. 
The time difference between the two instruments viewing the 
same Earth scene is within 30 s [14]. These orbits are selected 
over different latitude regions and seasons so that the collected 
data sets cover a wide range of scene temperatures. Sampled 
MODIS LIB 1-km and AVHRR LIB 4-km global area cover- 
age pixels within a 50-km area at nadir are averaged to match 
the horizontal resolution of the AIRS-retrieved atmospheric 
profiles (Fig. 2). To reduce the impact of spatial inhomogeneity 
on the average, sampled pixels within each 50-km area are 
examined for the standard error of T. Averaged points with 
standard errors larger than 2.0 K are excluded from our results. 

Because AVHRR channels have a relatively wider distri- 
bution of RSR than those of the spectrally matched MODIS 
bands (Fig. 1). it is expected that scene temperature and at- 
mospheric water vapor can have an impact on the differences 
between Xmodis and Tayhrr- Fig- 3 shows AT obs versus 
T for the 11.0- and 12.0-/zm bands. As expected, in general, 
AT is dependent on T. The 12.0-/zm band shows stronger 
dependence cn T at high temperatures. Recent studies used the 
split-window technique (i.e., the differences of T between 1 1 .0- 
and 12.0-/xm bands) to determine the total atmospheric water 
vapor content [15]. Thus, it is important to remove the impact 
of existing atmospheric water vapor on the calculation of AT 
in order to detect small calibration differences between sensors. 

Since the atmospheric temperature and water vapor profiles 
change constantly with time and geolocation, it becomes nec- 
essary to use simulation with real-time atmospheric profiles so 
that their effects on the calculation of AT can be captured. 
In this letter, the atmospheric profiles of retrieved temperature, 
water vapor, atmospheric pressure and ozone, and surface skin 


Fig. 3. AT obs versus T for the two spectrally matched bands/channels of 
Aqua MODIS and NOAA-18 AVHRR at 11.0 and 12.0 ^m. Each point is an 
averaged value over the AIRS pixel footprint at a nadir resolution of 50 km. 
Data are collected from Aqua MODIS and NOAA-18 observations collected 
from 13 SNO events in nearly overlapping orbits between 2006 and 2010. 



Fig. 4. AT simu versus T for the two spectrally matched bands/channels of 
Aqua MODIS and NOAA-18 AVHRR obtained from MODTRAN simulations 
with measured real-time atmospheric temperature, water vapor, pressure and 
ozone profiles, and skin temperatures from AIRS standard product. All data 
points come from simulations with atmospheric profiles matching in time and 
location with the SNO points shown in Fig. 3. 

temperature are taken from the AIRS standard retrieval product 
[16], collected at the same time over the MODIS and AVHRR 
SNO regions. The profile horizontal resolution is 50 km, which 
matches with the area in size used to average over for sampled 
pixels, and the vertical resolution is 28 levels between 1000 
and 0.1 mb. Given each AIRS-retrieved real-time atmospheric 
profile, a simulation of TOA radiance spectrum [L( A)] at a 
spectral interval of 0.001 /im (1.0 nm) is produced using 
MODTRAN4 [12]. The simulated spectral radiance L simu of 
the MODIS and AVHRR bands is determined based on (4). 
Then, the simulated Tmodis and ^avhrr are computed using 
(1) and (2). Fig. 4 shows A T simu versus T for the two MODIS 
and AVHRR bands at 11.0 and 12.0 /zm. Results show that there 
is almost no temperature dependence for AT simu when T is 
lower than 300 K. This is likely due to the fact that most SNO 
events for low-temperature scenes occur at around 70° N/S lat- 
itude, where the atmospheric water content is significantly low. 
Results also indicate that at high temperatures, AT simu shows a 
stronger dependence on T in the 12.0-/zm band than that in the 
11.0-/zm band. Comparison between observed and simulated 
AT provided in Figs. 3 and 4 reveals that the observed AT has 
a stronger dependence on T. 
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Fig. 5. Remaining temperature differences versus T obtained by subtracting the simulated AT from the observed AT. Three vertical lines represent temperatures 
at 0.3 of the typical radiance, the typical radiance, and 0.9 of the maximum radiance for MODIS bands. 



Fig. 6. Remaining differences in radiance (in percentage) versus T. 

B. Calibration Differences 

When a subtraction between the observed and simulated AT 
pairs of points is conducted, i.e., between Figs. 3 and 4, the 
remaining differences are considered to be caused by a relative 
calibration bias between MODIS and AVHRR according to (5). 
Fig. 5 shows the remaining differences versus T for the 1 1 .0- 
and 12.0-/um bands. To better understand the reason for the 
remaining differences related to calibration, the differences are 
examined at three temperatures corresponding to three radiance 
levels used in MODIS prelaunch calibration: 0.3 of the typical 
radiance, typical radiance, and 0.9 of the maximum radiance, 
respectively. Lines corresponding to these three temperatures 
are marked in Fig. 5. Values of the typical radiance are pro- 
vided based on the science product [1]. MODIS calibration 
requirements are 0.5% at the typical radiance for the 11.0- and 
12.0-/xm bands. At the 0.3 typical radiance and 0.9 maximum 
radiance, an extra 1.0% is added to the calibration requirements 
[17], indicating that the calibration quality is ensured for tem- 
peratures between T at 0.3L typ and T at 0.9 L max . For AVHRR, 
there is a similar stringent requirement at typical radiance, 
mainly focusing on the accuracy of sea surface temperature 
retrieval (http://www2.ncdc.noaa.gov/docs/klm/html/c7/sec7- 
l.htm). Results show that there is an excellent agreement be- 
tween Aqua MODIS and NOAA-18 AVHRR for the 11.0-/xm 
band for T from temperatures at 0.3L typ to 0.9L max , with a 
mean difference of less than 0.09 K. For the 12.0-//m band, 
the mean difference at L typ is 0.32 K, and there is a slight 
temperature-dependent trend with a warm bias (up to 0.60 K) 
for temperatures both close to T at 0.3T typ and 0.9L max . Fig. 6 



TABLE I 

Relative Calibration Differences and 1 -sigma Standard 
Errors Between MODIS and AVHRR 1 1.0- and 12.0-^m Bands at 
the 0.3 Typical, Typical, and 0.9 Maximum Radiances in Terms 
of Temperature (in Kelvins) and Radiance (in Percentage) 


Band 

0.3Ltyp (K) 

1.0 Ltyp (K) 

0.91. max (K) 

Average (K) 

1 1 .0 pm 

0.16 ±0.21 

-0.30 ±0.24 

-0.13 ±0.22 

0.09 ±0.42 

12.0 pm 

0.40 ±0.15 

0.36 ±0.21 

0.68 ±0.27 

0.32 ±0.42 


0.3Ltyp (%) 

1.0 Ltyp (%) 

0.9Lmax (%) 

Average (%) 

11,0 jim 

0.40 ±0.50 

-0.44 ±0.36 

-0.1 6 ±0.28 

0.39± 1.16 

12.0 pm 

0.93 ±0.38 

0.50 ±0.28 

0.81 ±0.31 

0.75 ±1.12 


shows the same differences shown in Fig. 5 in terms of radiance 
percentage, and results are also listed in Table I. Between T’s 
at 0.3L typ and 0.9L max , the differences are within 0.40% for 
the 1 1.0-/xm band, while for the 12.0-jxm band, the differences 
are 0.5% for T at L typ and up to 0.8%-0.9% for either T at 
0.3I/ typ or T at 0.9 T max . These differences are well within 
MODIS and AVHRR combined calibration requirements and 
indicate an excellent calibration consistency between the Aqua 
MODIS and NOAA-18 AVHRR bands at 11.0 and 12.0 /xm. 
It is also noticed that the remaining differences are temperature 
dependent, particularly for the 1 2.0-/xm band. A cold-scene bias 
is observed for temperatures lower than 240 K. At least half of 
the bias is due to MODIS calibration based on a comparison 
between MODIS and AIRS [18] and with in situ temperature 
measurements obtained in the Dome C area, Antarctica [19]. 
Improvements made in the upcoming MODIS Collection 6 LIB 
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product will significantly remove most of the cold-scene biases. 
For temperatures larger than 300 K, there might be land surface 
emissivity impacts if the emissivity has significant differences 
in spectral behaviors within the W-ll-fim range. From a 
calibration point of view, how well the nonlinear behavior is 
characterized has a direct impact on the calibration for low or 
high temperatures because the temperatures of on-board BB are 
limited to a small range. 

It should be pointed out that the differences found in this let- 
ter have minimum impacts on most skin temperature retrieval. 
This is because these retrieval algorithms are developed by 
regression against measurements such as buoy data [20], [21]. 
Thus, these algorithms are basically sensor dependent, and any 
constant calibration biases have been absorbed into the regres- 
sion coefficients. We also analyzed the sensitivity of simulated 
AT to errors in column atmospheric water vapor content. A 
15% change in water content only causes a AT change by less 
than 0.05 K at the typical temperature for a number of default 
atmospheric profiles provided by MODTRAN, indicating that 
impacts due to errors in retrieved water vapor product are 
generally negligible for cold to typical temperatures. 

IV. Conclusion 

This letter has compared the calibration differences between 
the Aqua MODIS and NOAA-18 AVHRR bands at 11.0 and 
12.0 fjbm using SNO observations obtained in nearly parallel or- 
bits. The temperature difference (AT) between the two sensors 
is used to evaluate their calibration at each scene temperature 
level from 0.3 of the typical radiance to 0.9 of the maximum ra- 
diance. The impact on AT due to the RSR differences between 
MODIS and AVHRR is estimated by MODTRAN simulations, 
and this effect is removed from observed AT. The simulations 
use real-time atmospheric profiles of temperature, water vapor, 
atmospheric pressure and ozone, and surface skin temperatures 
obtained from the Aqua AIRS standard retrieval product. The 
results of this letter show that the remaining temperature dif- 
ference is within 0.30 K (or 0.40% in radiance) across the 
entire temperature range for the 1 1 .0-//m band. For the 1 2.0-/jm 
band, the differences are 0.50% at the typical radiance and 
up to 0.90% at 0.9 of the maximum radiance. These results 
indicate an excellent calibration consistency between the Aqua 
MODIS and NOAA-18 AVHRR bands at 11.0 and 12.0 jUm, 
respectively. 
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